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The angle-resolved photoemission spectroscopy (ARPES) is able to measure both the spectra 
and spectral weights of the quasiparticles in solids. Although it is common to interpret the band 
structure obtained within the density-functional-theory based methods as quasiparticle spectra, 
these methods are not able to provide the changes in spectral weights of the electron excitations. 

We use Vienna Ab initio Simulation Package (VASP) for evaluation of the quasiparticle spectra and 
their spectral weights within Hedin’s GW approximation (GWA) for Fe^Si and a—FeSi 2 , providing, 
thus, a prediction for the ARPES experiments. Comparison of the GGA-to-DFT and GWA band 
structures shows that both theories reflect peculiarities of the crystal structures in similar shape 
of the bands in certain fc-directions, however, in general the difference is quite noticeable. We find 
that the GWA spectral weight of quasiparticles in these compounds deviates from unity everywhere 
and shows non-monotonic behavior in those parts of bands where the delocalized states contribute 
to their formation. Both methods lead to the same conclusion: those of iron ions in Fe^Si which 
occupy the positions Fe^\ where they are surrounded by only Fe ions ( Fe^ 2 ' 1 positions), have 
the d-electrons localized and large magnetic moment whereas Fe ions in the Fe 1 ' 2 ' 1 positions with 
the Si nearest neighbors have d-electrons delocalized and the magnetic moments quenched. The Si 
influence on the Fe ion state is even more pronounced in the a — FeSii where all iron ions have the 
Si ions as nearest neighbors: both GWA and GGA calculations produce zero moment on iron ions. 

The advantages and disadvantages of both approaches are discussed. 

PACS numbers: 71.20.Be, 71.20.Eh, 71.20.Gj, 75.20.Hr, 75.47.Np, 

71.20.-b, 75.10.Lp, 71.15.Dx, 71.70.Ch, 71.45.Gm 

I. INTRODUCTION 

A hope to use the electron spin additionally to the 
charge as an information carrier has led to a develop¬ 
ment of the spintronics. One feasible way to exploit 
the spin degrees of freedom is to synthesize such mag¬ 
netic semiconductors, which, on the one hand, should 
be magnetic at room temperature, and, on the another 
hand, should be easily integrated with existing semicon¬ 
ductor industry. Therefore, it is desirable that it should 
be Si basedi. A magnetic moment can be added by 
a transition-metal constituents. The other way is to 
use a magnetic metal for injecting of spin-polarized elec¬ 
trons into, say, Sz-based semiconductor. The technolo¬ 
gies which create magnetic epitaxial multi-layer films on 
Si, produce an interface, which contains the compounds 
of TSi , where T is a transition metal. This makes iron 
silicides compounds highly perspective materials both in 
bulk and film form and a detailed understanding of their 
physics is on demand. Recently the formation of sin¬ 
gle crystalline.FeaSj phases in the Fe/Si interface has 
been demonstrated by several groups^-A The theoretical 
understanding of the ground-state properties like cohe¬ 
sive and structural properties is achieved long ago via 
first-principle calculations based on the various realiza¬ 
tion of local-density approximation to density functional 
theor y 5 10 . 

However, the experiments, that measure differential 


information (not thermodynamics), require for their in¬ 
terpretation a knowledge of the single-particle- excitation 
properties. For example, all photoeffect-based measure¬ 
ments belong to this class. The photoemission spec¬ 
troscopy (PES) provides direct measurement of the en¬ 
ergy of electronic quasiparticles. Its extension, the angle- 
resolved PES (ARPES), which uses the synchrotron ra¬ 
diation facilities, allows for extracting also perpendicular- 
to surface momentum dependence of quasiparticle en¬ 
ergy. Further refinement, the laser-based ARPES pro¬ 
vides even better accuracy and resolution. Recently de¬ 
veloped method, the time-resolved two-photon photoe¬ 
mission (TR-2PPE) spectroscopyii, can monitor the 
state of an excited electron during the course of its trans¬ 
formation by laser-induced surface reaction. Probably, 
this is the only method capable to directly measure the 
quasi-particle’s life-time. The transport and tunneling 
experiments are even more evident examples where the 
qusiparticle concept is a necessary ingredient for under¬ 
standing the underlying physics. However, general theo¬ 
ries sometimes are not sufficient for describing real ma¬ 
terials. For example, the predictions of the life-times 
within the Landau theory of Fermi-liquid do not describe 
the experiments even on Al and noble metals^, contrary 
to the expectation. Indeed, electrons in these metals are 
well delocalized and expected to behave as a good Fermi- 
liquid. Furthermore, TR-2PPE experiments show that 
the lifetime of an excited electron in Al at a fixed energy 
E < Ep depends on the frequency of the pump pulse, i.e., 
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on the band from which the electron originated. These 
examples convince that ab initio calculations for real 
materials are required. The most developed method of 
electronic structure calculations, which is based on den¬ 
sity functional theory, is designed and applicable for the 
ground state properties onl y 12 ) 19 . This gives rise to cer¬ 
tain doubts that the electronic bands obtained within the 
framework density functional theory (DFT) by means of 
local-(spin)-density approximation with or without gra¬ 
dient corrections, can be interpreted as energies of ex¬ 
cited states. Nevertheless, the band structure, generated 
by Kohn-Sham equations, is ubiquitously exploited as 
quasiparticle spectrum.Here we will compare the results 
of calculations within the version of gradient-corrected 
local-density approximation to density functional theory 
(GGA) and GW approximation (GWA)—. In the latter 
G is the Green’s function and W is the screened Coulomb 
interaction. Both approaches have their own advantages 
and disadvantages. The exchange-correlation potential 
used in all modern implementaions of the Kohn-Sham 
machinery contains the correlation effects that are much 
beyond the random phase approximation used in GWA. 
Nevertheless the DFT based approach has the essential 
drawbacks that i) the exchange-correlation potential is 
calculated for the homogeneous electron gas, and the 
transfer of these results to the inhomogeneous electron 
gas of real materals is not easy to justify, and ii) it is dif¬ 
ficult to impove the calculations by adding in a controlled 
way some corrections. The convincing example of it is the 
widely used LDA-fU approximation, where many differ¬ 
ent forms of double-counting corrections are in use: the 
way to make it in a controlled way is not found yet. The 
GW method does not contain this problem, it is well con¬ 
trolled approximation. The strong advantage of GWA is 
that it results in real quasiparticle excitations. However, 
it requires much more of computer resources than, say, 
GGA. It is expected that they should produce approxi¬ 
mately the same band structure, since both approaches 
ignore the intraatomic (often strong) correlations. If the 
problem of interest does not require the knowledge of 
the spectral weights and/or is intended to use it as a 
first step for a better approximations, it is preferable 
to use GGA. The latter usually is used as a perturba- 
tional input for GWA. For this reason it makes sense to 
compare both quasiparticle and Kohn-Sham band struc¬ 
tures and to analyse the difference between the exchange- 
correlation potential and the real part of the GW self 
energy. Since the eigenvalues of the Kohn-Sham (KS) 
equations cannot be interpreted as quasi-particle (QP) 
excitations whereas the band structure generated by the 
GW method can, we calculate GW QP energies. They 
can serve as a predicton of ARPES measurements on iron 
silicides. 

The paper is organized as follows. In Sec. II we present 
the details of our ab initio DFT and GW calculations 
and give brief discussion of their advantages and draw¬ 
backs. The comparison of band structure and densitiy 


of electron states calculations of iron silicides Fe^Si and 
a — FeSii obtained by both methods is given in Sec. Ill 
and in Sec. IV we discuss results and make conclusions. 


II. THE METHODS: GW VERSUS GGA FOR 
IRON SILICIDES 

A. Band structure and excitation energies. 

The question if the DFT band structure can be in¬ 
terpreted as single-electron excitations’ energies have 
been discussed repeatedly (see, e.g.^, review^ and the 
bool—). The single-electron Green’s function by defi¬ 
nition describes the transitions between n and (nil) 
electron states, and the solutions of the Dyson’s equa¬ 
tion provide the excitation energies and their spectral 
weights. The potential in the Kohn-Sham equation does 
not depend on energy and automatically looses infor¬ 
mation on possible deviations of the spectral weight of 
quasi-particles from unity. Although the DFT based 
methods have been designed for the investigation of 
the ground-state properties and, it seems, there are no 
grounds for interpreting the band structure, generated by 
the Kohn-Sham equations as the energies of the single¬ 
electron excitations. However, the question about such 
a possibility is not so trivial and cannot be ruled out 
by this argument only. First, one can notice that the 
same statement holds for the Hartree-Fock-determinant 
wave function of the ground state , for which Koopmans^S 
proved that in closed-shell Hartree - Fock theory, the 
first ionization energy of a molecular system is equal 
to the negative of the orbital energy of the highest oc¬ 
cupied molecular orbital. There is an analogue of the 
Koopmans’ theorem for DFT—. Second, Hoenberg-and- 
Kohn’s proo^. that the total energy of a system with 
interaction is a unique functional of the electron den¬ 
sity p( r) led Sham and Schlueter— to the equation, 
which, at least, in principle, allows calculation of the 
exchange-correlation potential in DFT for any chosen ap¬ 
proximation for the electron-Green’s-function self-energy 
S'. Schematically, G _1 = (Gq 1 — v xc ) + v xc — E' gives 
G = Gdft + Gdft ( v xc — E') G. Here Go is the Green’s 
function of electrons in the Coulomb field of nuclei and 
Hartree potential and the prime in E' means that the 
Hartree term is subtracted, E' = E — vHartree- Express¬ 
ing then the charge densities in terms of G, p = ( G ) E , 
i.e., taking appropriate contour integral over energy from 
both sides and using the equality p^J d {r) = p ge nuine{r) 
one finds 0 = ( Gdft (v X c — E') G) E . This equation^ es¬ 
tablishes the correspondence between the self-energy and 
the exchange-correlation potential. For this reason one 
can hope that at least for the v xc found for certain E' 
the self-energy-based and Kohn-Sham eigenvalues coin¬ 
cide. The problem is, however, that the eigenvalues £k n 
and wave functions (r) of Kohn-Sham equation, 
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represent more differential information than the elec¬ 
tron densities and, to the best of our knowledge, a 
prove of a statement that e kn = E kn does not ex¬ 
ists ( here n labels the bands). To be more pre¬ 
cise, Eq. m is an approximate version of the equa¬ 
tion for the functions $fc n (r, w), which diagonalize 
the Green’s function in GWA at arbitrary energy 
( < hfcn(^ J )| G (ca) [^^^(cc)) = d kk '6 nn > A kn (ca)), 

i.e., it contains energy-dependent self-energy E'(r, ri,w) 
and $fc ra (ri,w); the eigenvalues, which are poles of the 
Green’s function, are obtained from the equation E kn = 
ReA kn ( E kn ). At last, both equations are highly non¬ 
linear and uniqueness theorem does not exist for them. 
Particularly, the problem of choice of the physical so¬ 
lution from multiple solutions generated by GWA is dis¬ 
cussed in ref4i. Thus, at the moment we are at the stage 
when even the form of the exchange-correlation potential 
for a self-energy in GW approximation is not derived and 
the question about relationship between Kohn-Sham and 
GWA energies remains open. It is interesting, however, 
that more accurate exchange-correlation functionals than 
the widely-used ones can be obtained for evaluation of 
the excitation energies: the authors of ref£ state that, 
on the one hand, a surprising degree of agreement be¬ 
tween the exact ground-state Kohn-Sham eigenvalue dif¬ 
ferences and excitation energies, for excitations from the 
highest occupied orbital to the unoccupied orbitals has 
been achieved for small systems; on the other hand, that 
“ the popular LDA and GGA functionals are nowhere near 
sufficiently accurate”—. 


B. Band structure and spectral weights. 


the excitations. Contrary to that the self-energy does 
depend on energy giving birth to the spectral weight of 
the GW-QPs Z e < 1 ; the remaining part of the spectral 
weight is shifted to incoherent excitations. It is instruc¬ 
tive to represent the spectral weight in terms of magni¬ 
tudes, encoded in the VASP. Let us write the equations 
for the Green’s functions, corresponding to Eqs. ©,©: 

G~\E) = G^ ft (E) - (E '(E) - v xc ). (3) 

Within the “one-shot” (GoWo) approximation E = Ek n 
and we obtain 
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This relation is fulfilled by the magnitudes 
Ekni £km (vxc ) , E Gff (e b ), which are gener- 
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ated by VASP. As seen from Eq.© in the fori 
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the role of the coefficient in the first term of the 
quasiparticle-energy expansion with respect to presum¬ 
ably small perturbation ReT, ( e kn ) — v xc , is played by 
the spectral weight. If the perturbation theory works, 
E kn - £kn < (v xc - E) fen everywhere. Our calcula¬ 
tions confirm that this inequality is valid for iron sili- 
cides. A comparison of behavior of these two differ¬ 
ences, fc-dependence of the spectral weights and an anal¬ 
ysis of contributions from different electrons into these 
fc-dependences are given in Sec. dHI CD . 


The PES-based experiments provide also the informa¬ 
tion on the spectral-weights of QPs, Z E . As was dis¬ 
cussed above, the “self-energy” in Kohn-Sham equation 
m does not depend on energy and, therefore, corre¬ 
sponding spectral weight, defined in the perturbation 
theory as (Z k \ = [1 — dT,(E)/dE} -1 ^ , is always 

equal to unity. Insufficiency of the DFT-based approach 
is seen also from the Sliiter & Sham’s equation for the 
exchange-correlation potential: since it is based on the 
equalities of model and genuine electron densities of 
many-electron system, which are integral characteristics, 
it does not provide the equation for spectral weights of 


III. IRON SILICIDES 

The calculations presented in this paper are per¬ 
formed using the Vienna Ab initio Simulation Pack¬ 
age (VASP)^ with Projector Augmented Wave (PAW) 
pseudopotentialsdd. The valence electron configurations 
3d 6 4s 2 are taken for Fe atoms and 3s 2 3p 2 for Si atoms. 
One part of calculations is based on the density func¬ 
tional theory where the exchange-correlation functional 
is chosen within the Perdew-Burke-Ernzerhoff (PBE) 
parametrizationi^ and the generalized gradient approxi¬ 
mation (GGA). 
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In the GW part of calculations implemented in VASP— 
we report only one-shot approximation, so-called Go Wo- 
The GGA Kohn-Sham band structure and eigenfunc¬ 
tions were taken as the input for the GW calculations. 
The self-energy £ is computed as £ ~ iG GGA W GGA . 
Throughout both GGA and GWA calculations the plane 
wave cutoff energy is 500 eV, the Brillouin-zone integra¬ 
tions are performed on the grid Monkhorst-PacI— special 
points 10 x 10 x 10 for Fe^Si and 12x12x6 for a — FeSi2 
and Gauss broadering with smearing 0.05eV is used. For 
all GW calculations the number of frequencies is 500 and 
160 electronic bands are used. We use the complex shift 
of 0.043 eV for the calculation of self-energy. 

A. Iron Silicides structure 

The calculations are performed for two iron silicides, 
Fe^Si and a — FeSii- The compound Fe^Si belongs 
to DO 3 structural type with the space symmetry group 
Fm2>m. The iron atoms have two nonequivalent crys¬ 
tallographic positions in fee lattice, namely, Fe^and 
Fe have different nearest surroundings: Fe W has 
eight Fe^ nearest neighbors which form a cube, whereas 
the Fe^ is in the tetrahedral surrounding of both Si and 
Fe^ atoms. 

The iron bi-silicides have several structural modifi¬ 
cations. The most stable phases are a — FeSi 2 and 
ft — FeS'G phase s 23 ) 24 . The compound a — FeSi 2 has 
tetragonal lattice with Pd/mmn space symmetry group 
with one molecula per unit cell. Each iron atom here 
is located in the center of cube, consisting of the sili¬ 
con atoms. This structure contains the planes which are 
formed by only iron and by only silicon atoms. These 
planes are orthogonal to tetragonal axis. Two planes 
formed by silicon atoms are separated by wide empty 
cavity, which does not contain the iron atoms. Accord¬ 
ing to our calculations this cavity is expected to play 
essential role in the transport properties of a — FeSi 2 - 

The rhombohedral cell have been used for the Fe^Si 
calculations. The equilibrium parameters and the dis¬ 
tances between nearest Fe and Si atoms for the Fe^Si 
and a — FeSi2 structures have been found from the full 
optimization of the structure geometries within GGA and 
are shown in Tab. ©■ The GW calculation have been 
performed with the same structural parameters. 

Both spin-polarized GGA and GW result in metallic 
states, ferromagnetic for Fe^Si and paramagnetic with 
zero-spin Fe atoms for a — FeSi2- For this reason all 
further calculations for a — FeSi 2 have been performed 
within non-spin-polarized version of VASP. The struc¬ 
tural inequivalence of the Fe-atoms’ surroundings in the 
Fe^Si reflects itself in both magnetic moment values and 
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As was mentioned above, different chemical surroundings 
of the Fe atom positions, the cubic one for Fe W by Fe*- 2 -* 


the contributions of Fe-ions’ d-states into electron den¬ 
sity of states. The magnetic moment Mp e i of FeW atom 
is higher than the free-atom moment, M GGA = 2.52/is 
and M G ^ A = 2.55 ps what corresponds approximately 
to the d 5 configuration. The Fe^ atom has much lo wer 


Fe 3 Si 

a — FeSii 

a = 5.6A (5.65A) 

R (Fe (1) - Fe (2) ) = 2.45A 
R (Fe (2) - Si) = 2.45A 

R ! Fe (1) - Si) = 2.83A 

a = 2.70A (2.69A) 
c = 5.13A (5.134A) 
zsi = 0.272 (0.28) 
R(Fe- Si) = 2.30A 
R (Si — Si) =2.56A 


Table I: Relaxed lattice parameters and the equilibrium dis¬ 
tances between nearest ions. The experimental values^ are 
given in brackets. 


moment, M GG 2 A = 1.34 pb and M G ^ A = lAOps- As 
will be seen from the analysis of DOS, the latter mo¬ 
ments are formed by the delocalized d-states. The ex¬ 
perimental values reported in works^ and2& are slightly 
different: = 1.2 ps, = 2. dps in Ref— and 

Mpf 2 = 1.35 p B , MpJl = 2.2 /is in Ref.-^. 


B. Comparison of GGA and GW densities of 
electron states 

Fig. dT]) displays comparison of the GGA and GW den¬ 
sities of electron states (DOS) for Fe^Si and a — FeSi 2 . 
The GGA part of results coincides with previous calcu¬ 
lations of FesS a 25 , 27 , 28 and a — FeSi^r^r^. The gen¬ 
eral features of the DOS in both compounds and approx¬ 
imations are that the bands in the interval (—5, +5) eV 
around Fermi energy are formed by the d-electrons of iron 
with a slight admixture of s— and p— electrons of Si and 
Fe. The Si valent s- and p —electrons are delocalized in 
the wide energy region with smeared maximum around 
—4 eV in both compounds. GWA changes the intensities 
of the peaks mainly in the energy region deeply under 
Fermi surface, but the changes in Fe^Si and a — FeS%2 
are different. If in the GGA DOS of the Fe^Si the peak 
located at F ~ —3.5 eV is shifted by GWA for about 0.5 
eV and made sharper, the GGA peaks in a — FeSi2 DOS 
in approximately the same energy region (I, II in bottom 
panel of Fig. (JT])) is washed out within the GW calcula¬ 
tions. There is a difference in the GWA vs GGA changes 
in DOS for spin “up”(majority) and “down”: while the 
GW spin-majority states remain almost untouched com¬ 
pared to GGA, the spin-minority GWA peaks are shifted 
towards the Fermi energy. 


atoms and the tetrahedral one for the Fe^ atoms by the 
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Figure 1: The left panel: The spin-polarised density of electron states for Fe^Si in the interval of energies [-4eV, 2eV]. Since 
the GGA and the GW approximations produce different Fermi energies, £ GGA (FesSi) = 7.88 eV and e GW (Fe^Si) = 8.44 eV, 
the plots are aligned for comparison by placing the zero in the energy axis of both plots at Fermi energy. The right panel: The 
DOS for a — FeSi 2 in the energy interval [-5,5]eV with the same type of alignment of the energy axes; £ GGA (a — FeSi 2 ) = 9.34 
eV and £$ w (a - FeSi 2 ) = 10.03 eV. 


Fe 1 ' 1 ' 1 and Si atoms as nearest neighbors reflect them- is illustrated in Fig. ([21). 
selves in different behavior of partial d-electron DOS. It 

I 


The general shape of partial d-DOS is not changed 
by GWA. The d-states of E g symmetry of Fe ' 1 ' 1 as well 
as Fe do not contribute into DOS at the Fermi level 
in both in LDA and GWA (panel (A), (C), Fig.(0). As 
seen from Fig. ([2]), panels (A),(B), the contribution to the 
magnetic moment on Fe^ from E g -orbitals (positive 
DOS in (A)) compensated by the contribution from T 2g - 
orbitals ( the negative middle peak in (B)). This means 
that the T 2g -orbitals are responsible for formation of the 
large quasi-localized magnetic moment at Fe' 1 ' atoms. 
It is also interesting that the usual splitting of the d-shell 
into T 2g - and E g - symmetries is violated here by the con¬ 
tribution from exchange interaction: as seen from (A) 
and (B) panels the F g -peak is in-between two T 2g -peaks. 
A different behavior is demonstrated by d-electrons of 
Fe™-. its d-DOS is spead in a wide region of energies. 
The d-electrons of both T 2g - and E g - symmetries con¬ 
tribute to formation of magnetic moment of Fe^ 2 \ The 
delocalization of Fe^ d-electrons reflects itself in the 
smaller moment than the one on Fe W. The d-DOS of 
iron in a — FeSi 2 , where Fe atoms also have Si atoms as 
neighbors, displays behavior similar to Fe ' 2 ' 1 partial DOS 
for d-electrons in Fe^Si . However only these delocalized 
electrons are not able to form magnetism in a — FeSi 2 . 
Indeed, the criterium for magnetism formation for the 
delocalized electrons is much harder to fulfil than for the 
case of localized electrons whose role is played by T 2g - 
electrons of Fe 1 ' 1 ' 1 in Fe^Si . The absence of magnetism 
in a — FeSi 2 is easy to understand on the basis of well- 
known Stoner’s model for a magnetism of the delocalized 
electrons: the criterium Jg{ep) > 1 is not fulfilled since 
the density of electron states <7 (ef) at the Fermi energy 
£f is too small (here J is exchange integral between delo¬ 
calized electrons). An alternative mechanism of the mag¬ 


netism supression would be a formation of the low-spin 
state within the localized d-electron picture. This state 
could be formed if the crystal-field splitting of the d-shell 
was stronger than the Hund-exchange one. However, the 
density of d-electron states does not contain bright peaks 
which might be interpreted as former d-levels splitted in 
the crystal field. Thus, one can conclude that if GGA 
and GWA are good approximations for a — FeSi 2 , the 
key mechanism of the magnetism destruction in this com¬ 
pound is the delocalization of d-electrons. 

The most pronounced changes in GWA compared to 
GGA are experienced by T 2 9 electrons. It is illustrated 
on Fig. ([2]) , panels (E,F) for a — FeSi 2 : two peaks (I and 
II) seen in the GGA DOS which have been mentioned 
above and washed out in the GWA are formed namely by 
T 2g electrons. The same is valid for the “down” spin T 2g 
states in the vicinity of the Fermi level in Fe^Si (panel 
(D) in Fig. (0). At the same time the well-expressed lo¬ 
calized peaks formed by E g orbitals remain intact. This 
is understandable: the GWA is the approximation which 
takes into account the effects of screening which is pro¬ 
vided namely by the delocaized electrons. These effects 
are taken into account within the GGA and GWA differ¬ 
ently and, therefore, one can expect that the difference 
will be more pronounced for these types of electrons. 


C. Comparison of GGA and GW band structures 
and analysis of the spectral weights 

The Kohn-Sham band structure calculated within 
GGA does not differ from the known results for 
Fe^S i 2828 and for a — FeSi^,— ■ Here we com¬ 
pare the GGA and GW bands in some of symmet- 
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Figure 2: The partial spin-polarised d-electron DOS for FezSi (upper four panels, (A),(B),(C),(D)) and a — FeS^ (lower two 
panels, (E),(F)) . Left panels display the contribution to DOS from d z 2 — and d x 2_ y 2 — states ( E g ), while the right ones show 
the contribution from d xy — , d xz — and d yz — states (Thg). 


ric directions. Figs.© and © show the band struc¬ 
ture for Fe^Si in the directions TX and TL, and for 
a — FeSi 2 in the directions TX , TM and TZ 1 where 
r = (0,0,0), X = (27T/a) (1,0,0), L = (tt/o) (1,1,1), 
M = (2-k/o) (1,1,0), Z = (27r/c) (0,0,1). 

We report here the results of the comparison only for 
the part of the GGA and GWA band structures which are 
within several electron-volts vicinity of the Fermi energies 
(remind that GGA and GWA generate different Fermi en¬ 
ergies, see the capture to Fig. ©)• The bands are named 
in accordance with their symmetries in the T—point. The 
closest to the Fermi energy three filled spin-up bands for 


the Fe^Si in Fig.© left) near T—point, the dublet E g 
and the triplet T 2 gi are formed by the d-electrons of Fe 
atoms. First empty band ( A\ g ) near F—point is formed 
by the s-states of both Fe and Si atoms. 

The GGA and GWA band structures for a — FeSi 2 
are shown in Fig.©. Here the closest to the Fermi 
energy filled bands are formed by the d-orbitals of Fe 
atoms near T— point (T4/1 group) are Bi g ( d x 2-^2), A\ g 
(d z 2), the dublet E g ( d xz ,d yz ) and i?2g( d xy ). The low¬ 
est shown band (A. 2 „) is formed by the s-electrons of Fe 
and p-electrons of Si. Again we observe the same ten¬ 
dency: namely the delocalized states, in this case, s- and 
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Figure 3: The GGA and GW spin-up (top) and spin-down 
(bottom) bands for FezSi. Zero in the energy axis of GGA 
and GWA plots is chosen at corresponding Fermi energies. 


p-states, show the largest difference in GGA and GWA. If 
the band formed by d-electrons close to the Fermi energy 
remain almost untouched, the lowest sp band is shifted 
in GWA by ~ 1 eV. First empty band (E u ) near F—point 
is formed by the p-states of both Fe and Si atoms. In 
general the GGA vs GWA shift is around half of electron- 
volt for the excited states, while the band shape remains 
the same. As seen at the right panel of Fig. d4j in the TZ 
direction, the purely d-bands are completely flat, while 
the dispersion which arises near the boundaries is due to 
the admixture of s and p states. Analogous admixture 
of sp electrons is observed around the boundary points 
X and M. This is easy to understand since the motion 
along TX direction in the £:-space corresponds to motion 
from the Fe atoms plane to the plane of Si atoms in real 
space; Si atoms do not have d-electrons and can admix 
the s- and p-states only. 

Here we return to the question asked in Sec. (Ill Bl) . 
which electrons mainly are responsible for the deviations 
of GWA from GGA results. Let us consider the example 
of two Eg and T 2 g filled bands of quasiparticles for spin 



Figure 4: The band structure of a — FeSi 2 ', zero in the energy 
axis of GGA and GWA plots is chosen at corresponding Fermi 
energies. 


“up” shown in Fig.©. 

First general feature seen from both Fig.® and Fig.© 
is that the spectral weights Z kn , that within the Kohn- 
Sham approach are equal to one by construction, within 
the GWA are strongly decreased. Second feature, some¬ 
what surprising, is that the non-monotonous behaviour 
of Z^n arises everywhere where an admixture of delocal¬ 
ized electron states is present. Indeed, one observes that 
near the X point of the Brillouin band it is the admix¬ 
ture of s- and p-electrons to d-states for the band T,^ 1 
of Fe 3 Si causes the changes in the Z kn ; the difference 
ReYi' — v xc is changed faster than E kn — £kn (see panel A 
at Fig.©). The picture for the empty bands is different: 
the s- and p-states of Si and d-states of Fe are mixed 
in the center of the band, whereas the contribution of 
the d-states is increased around boundary points X and 
L. The quasiparticle energies of the excited states are 
lower than their Kohn-Sham counterparts. Again, the 
largest difference is observed in the those parts of the en¬ 
ergy spectrum where the contributions from s— and p— 
states becomes significant. 

As seen from Fig.©, the spectral weights Z kn oi the 
quasiparticles in a — FeSi 2 also show a strongly non- 
monotonous dependence on k. For example, the spectral 
weight Z k ,B lg of the Bi g band in the direction TX (sec¬ 
ond from top panel in Fig.©B) shows sharp changes in 
the interval 0.55 < Z k ,B lg < 0.8 due to much faster de¬ 
pendence on k of the difference (S' — v xc ) than of the 
( Eqp — egga ) one. The lower part of the panel B in 
Fig.© explains the reason: again, the closer to the X 
point the higher contribution from the delocalized s- 
electrons. This confirms, thus, the general tendency no¬ 
ticed above. 


IV. DISCUSSION AND CONCLUSIONS 

The comparison of the band structures obtained in 
the ab initio calculations within the VASP for Fe^Si 
and a — FeSi 2 in GGA and GqWqA shows that in gen- 
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Figure 5: Fe^Si quasiparticle filled bands (top panels), their 
spectral weights Zk (second from the top), the fc-dispersion 
of numerator and denominator of Eq.® (third from the top) 
and the character coefficients Ck (bottom panels, see text). 
The letters label different characters: d and p stand for the 
empty and filled circles correspondingly, the black squares 
denote s character. The filled bands are shown in (A) and 
(B) panels, while (C),(D) display the empty ones. The dashed 
and solid lines on (A) panel denote the non-degenerate bands 
Eg. The dashed and solid lines on (B), (D) panel denote 
non-degenerate and double degenerate bands Tbg- 


eral the bands shape is similar. The difference between 
GGA and GWA bands becomes more pronounced in 
those parts of the Brillouin zone where the delocalized 
states give noticible contribution into quasi-particle en¬ 
ergies. This observation is somewhat unexpectablc since 
both approximations are designed for description of well 
delocalized (Fermi-liquid like) electrons. There are at 
least two sources which can contribute to this difference. 
First is the fact that the standard GGA is not free from 
the self-interaction while GWA takes into account Fermi 
statistics by construction. Second source is that GGA 
and GWA are quite different approximations, as was dis¬ 
cussed in Introduction. 



Figure 6: The a — FeSi 2 bands. The notations are the same 
as in Fig. ©, but the filled bands are imaged at the panels 
A,B,C, while the empty one is shown at D panel. 


The other astonishing moment following from GWA cal¬ 
culations is that in spite of the fact that the electrons 
in both systems are well delocalized the spectral weights 
Zkn of the quasiparticle bands Ek n are decreased almost 
by a half. One could assume that the reason for that is 
that the basis set is not sufficiently large, however, the 
results are not changed with further increase of number 
of bands which are taken into account. This means that 
the GWA indeed shifts the remaining part of the weight 
to the incoherent part of excitations. 

Both GGA and GWA band structures and, corre¬ 
spondingly, the density of electron states, show that the 
d-electrons of those Fe atoms which have Si nearest 
neighbors, namely, Fe^ atoms for Fe^Si and all 
Fe atoms in a — FeSi 2 , are more delocalized than 
the d-electrons of Fe^ atoms in Fe^Si which have 
only the other Fe atoms as neighbors. The partial 
density of states of Fe W d-electrons with E g and T 2 g 
symmetry in the T point has well-expressed peaks, 
the positions of which could be ascribed to a splitting 
in the crystal field. However, this splitting does not 
correspond to the standard picture of the quasiatomic 
levels a d\ a d t a + e° e a d\ a d ea: from which the bands 
are formed, the interactions renormalize these “levels” 
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£ T 2g a £ t 2 g ,<T in such a way that their sequence 

becomes £t 29 ,-|- < £ t 29 ,4.j £e 9 ,t < £ T 2g ,f. 

We performed the calculations within the one-shot 
GW approximation (GoWo )> and it is worth to notice 
why this approximation is preferable compared to the 
fully self-consistent GWA. First, it often gives the 
results which describe PES experiments better (see, 
e.g., refs i 32 ) 33 ) 34 ). Second argument is connected with 
the nature of the self-consistency. Each term in the 
self-consistent perturbation theory (scPT) corresponds 
to whole series in the non-scPT. If we consider these 
corrections, we find that a part of corrections in each 
order of the non-scPT changes the matrix elements of 
the Coulomb interaction, actually, via corrections to the 
wave functions. The self-consistency loop changes the 
eigenfunctions which diagonalize the GW equation for 
the Green’s functions, and, correspondingly, dresses the 
matrix elements of Coulomb interaction by the RPA 
(random phase approximation) graphs only. All other 
contributions to the wave functions, which arise in the 
same orders of non-scPT, and, correspondingly, to the 
matrix elements of Coulomb interaction, are not taken 
into account. However, in the case of non-homogeneous 
electron gas in real-material, where the gellium-like parts 
can coexist with strongly correlated liquid, the situation 
it not alway obvious. The question how and to which 


terms the vertex corrections should be applied remains 
highly non-trivial. This have been demonstrated, e.g., 
in ref. 3 -, where it was shown that the vertex corrections 
into effective interaction W may improve the results 
while turning on the corrections in the self-energy may 
lead to unphysical quaiparticle dispersion. At last, the 
calculations within GWA with the vertex corrections 
demand for much stronger computer resources than 
even scGWA. The ARPES experiments on iron silicides 
would be of great help in further understanding of these 
compounds and, possibly, could motivate more advanced 
theoretical approach to the problem. 
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